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Abstract 

Avalanching systems are treated analytically using the renormalization group (in the self-organized- 
criticality regime) or mean-field approximation, respectively. The latter describes the state in terms of the 
mean number of active and passive sites, without addressing the inhomogeneity in their distribution. This 
paper goes one step further by proposing a kinetic description of avalanching systems making use of the 
distribution function for clusters of active sites. We illustrate application of the kinetic formalism to a model 
proposed for the description of the avalanching processes in the reconnecting current sheet of the Earth 
magnetosphere. 
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I. INTRODUCTION 



Many natural systems, which work in an open configuration, respond to external disturbances 
showing scale-invariant discrete events Jl|. One common feature of these systems is the develop- 
ment of a local threshold instability in an avalanching manner. In the late 80s the concept of self- 
organized criticality (SOC) was proposed by Bak et al. [2] for the dynamical- statistical behavior of 
such systems. SOC has been applied to a variety of systems (see Jensen |Q] and references therein 
for a list of some such systems). Although SOC can exist, strictly speaking, only in the limit of 
infinitely slow external input where complete separation of time scales is achieved [4] , it has also 
been applied to presumably avalanching systems with strong driving. A good example of such sys- 
tems is space plasma and, in particular, the plasma in Earth's magnetotail under magnetic substorm 
conditions |Q]. Since SOC is questionable for such strongly driven systems we, in what follows, 
address to them as to avalanching systems, bearing in mind the avalanche-like propagation of local 
instabilities. Up to date the most often used tool for studies of such systems is numerical modeling. 
The usual analytical approaches proposed so far are the renormalization group methods (see,e.g., 

nn Jill 

Refs. ylla] and references therein), and the mean-field description (see,e.g., Refs. y|,|7D and ref- 
erences therein). The renormalization group methods assume scaling from the very beginning and 
are applied only in the close vicinity of the stationary (critical) point, that is, in the self-organized 
criticality regime. The mean-field approach is based on the analysis of the mean number of active, 
passive, and critical sites. It is not restricted to the criticality range only including it as the limit 
of zero number of active sites. Mean-field approximations predict self-organized criticality in the 
limit of zero average number of active sites and, strictly speaking, are applicable only for system 
dimension exceeding some critical number, often well above the dimension of real physical sys- 
tems [8]. Mean-field obtained exponents are often consistent with those found experimentally and 
numerically for lower dimensions too but no quantitative explanation is given. On the other hand, 
deviations from these exponents for real systems are quite usual. The mean-field approach does 
not take into account the tendency of the active sites to organize in clusters. Indeed, if avalanches 
of various durations and sizes are present, the distribution of active sites at any moment should be 
very inhomogeneous. In the present paper we propose a novel approach to the analytical descrip- 
tion of avalanching systems which is based on the kinetic equation for the distribution function 
for active site clusters. We demonstrate the power of the kinetic formalism applying to the model 
which was recently proposed as a model of avalanching reconnection in the current sheet of the 
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Earth magneto sphere 



H. KINETIC EQUATIONS FOR CLUSTERS 

The mean-field approach has the obvious drawback of ignoring that active sites have the ten- 
dency to appear in clusters. These clusters are, in fact, the instantaneous snapshot of the developing 
avalanches, so that the size of each cluster is time dependent, w = w(t). However, when consid- 
ering many coexisting clusters, we may describe their behavior with the help of the distribution 
function f(w,t) = dN/dw, where now the cluster size w and time t are independent variables. 
The evolution of the single cluster size will be translated into the evolution of the distribution 
function. The total number of active sites is given by the integral 

/■oo 

N a = / wf(w)dw (1) 

J 

We have to introduce also the number of passive sites N(0) (similar to what is done in case of 
a Bose-gas, where the number of particles in the lowest state is macroscopically large). Then 
N a + JV(0) = const. 

Let P + (wi, w 2 ) be the probability of the cluster growth (per unit time), and P_(w 1 , w 2 ) be the 
probability of shrinking. Then 

df(w) 
dt 
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P+{w,w')f{w')dw' + / P-(w,w')f(w')dw' 

w (2) 
P+(w',w)f(w)dw' - [ P-(w',w)f(w)dw' + j(w)N(0) - ^ 



'o • • r(w) 

The term r y(w)N(0) in © describes the birth of active states due to external driving, while the 
last term takes into account the finite life-time of clusters, i.e. the transition to the passive state 
(Bose-Einstein condensation). If the driving is sufficiently strong and avalanche merging is not 
negligible, the kinetic equation should be completed with the with the time-dependent "non- 
linear" merging terms 

'df(w) 



dt 



Pi(w, Wi, w 2 )f(w 1 )f(w 2 )S(w — Wi — w 2 )dwidw 2 

(3) 

Pi(wi, w, w 2 )f(w)f(w 2 )S(w 1 — w — w 2 )dwidw 2 . 



Merging becomes progressively more important when the average fractional density of active sites 
increases. When this density is not too large (it does not have to be small though, in contrast with 



the SOC regime), merging will be still relatively weak and can be further studied perturbatively. 
Strong merging corresponds to the very strong driving, so that the system behavior is, at least 
partially, forced externally. In the present paper we assume that driving is moderate (not weak and 
not exceptionally strong) so that merging can be ignored at this stage, deferring treatment of very 
strongly driven systems to elsewhere. In our case one can expect that there is a wide range (inertial 
interval) of cluster sizes in which the distribution shape is independent of the external driving and 
is determined by internal dynamics and/or space dimension. 

In general, the distribution function f(w) would depend on the growth and shrinking probabil- 
ities. We shall consider here the class of systems where growth and shrinking occur only at the 
boundaries of clusters. It should be noted that the dynamics inside clusters may induce transitions 
between active and passive sites, producing, e.g. "punctuated" clusters for the classical sandpile 
model where an active site becomes passive at the next step. We shall measure the size of 
such cluster including the passive (receiving) sites as well, so that the internal dynamics does not 
affect the cluster size. Situation may be more complicated when clusters are developed fractals, 
with tunnels appearing and crossing the cluster JjJJ _ Such systems would probably require special 
treatment. We restrict ourselves here with the clusters which grow of shrink at their boundaries. 
Space and laboratory plasma systems liol seem to belong to this class. 

In this case the probabilities are nonzero only for \w' — w\ = A -C w, so that © can be written 

df 

— = -P_(w)a(w)f(w) - P + {w)a{w)f{w) 

dt (4) 
+ P_(w + A)a(w + A)/0 + A) + P + (w - A)a{w - A)f(w - A), 

where a(w) is the density of states. This approximation is not valid for small w, where the cluster 
kinetics should be strongly affected directly by driving. We seek for an approximate description of 
the cluster kinetics in the range where it is determined but the internal features of the system rather 
than by external influence. It is obvious, that if a large size strong driving is applied the reaction 
of the system would be a forced reaction and not self-organized in any way. 

The approximation may be not accurate for largest clusters either, since possible fractality 
111 of clusters may result in the breakdown of independence of probabilities at neighboring active 



boundary sites. Indeed, all numerical simulations [3] show distortions for very small and very large 
w. Thus, the physical sense of our approximation is that we are working in the inertial interval far 
from both limits. According to existing analyses, such interval exists almost always. 
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For one-dimensional clusters a(w) = 1 or a(w) = 2 (the latter holds for growth in both 
directions). This allows immediate n-dimensional generalization. Let w be a linear measure of 
a cluster (effective radius), and let D be the cluster volume. The density of states a is then the 
cluster surface area. In general, D oc w^, a oc vf , n>fj l >u>n — 1, where ji and v are fractal 
dimensions of the cluster volume and boundary, respectively. Taylor expanding © we arrive at 
the following differential equation 

where a = A(P_ - P + ) and f3 = (A 2 /2)(P_ + P + ). The stationary solution 



/ oc (l/pa)exp 



[a//3)dw 



(6) 



exists only if a > 0. In general, a and (3 can depend on w. Both describe the local growth (shrink- 
ing) per site at the cluster surface. Their dependence on the cluster size would mean essentially 
that the growth and shrinking probabilities as well as variation of the affected neighbor zone at 
some site depend on what happens at other sites. While, in principle, this cannot be excluded 
(waves could transfer information across the cluster or long range forces are involved yj]), many 
avalanche systems seem to be governed by local dynamics, so that it is natural to consider (at least 
at this stage) the case of probabilities independent of w (see, however, comment in section ITTTb. 
One finds 

/ oc o~ x exp(— w/w c ) = w~ v exp(— w/w c ), (7) 

with w c = const. The obtained / = dN/dw describes the distribution of linear sizes (effective 
radii). For the distribution of the cluster volumes one has 



dN dw dN 



oc ^a-M-")/!* exp f-AD^ . (8) 



dD dD dw 

In the mean-field limit n > 1 |Q] one has (dN/dD) oc D~ 2 . 

The derived expressions assume isotropy. If the system is anisotropic and/or a preferential shape 
of cluster exists, e.g. clusters are elongated Jl], the above treatment may have to be modified 
by considering vector w describing linear sizes along principal axes. These modifications are of 
technical character and do not change substantially the basic equations and conclusions. Yet, they 
require a more lengthy analysis and cannot be presented in a letter. We will provide this analysis 
elsewhere. 
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III. BURNING MODEL 

The above theory can be illustrated on the simple "burning" model J^] described below. In this 
model each site is characterized by its temperature, T(x). The external driving is random heating 
of the sites. The amount q of heat per unit time is going to a site with probability p, so that the 
average heat transfer from outside (in driving) is pq. The temperature of a passive site (the one 
which is not burning) changes according to 



where 77 (t) is a random number, \rj\ < 1, so that (rj(t)rj(t')) = 5(t — t'). Once T > T c , where T c 
is some critical temperature, the site becomes active. An active site burns and produces heat at the 
rate J = fxT, /i < 1. During the burning stage the temperature decreases (unless driving is strong 
enough to force permanent burning). When the temperature drops below some value TJ such that 
T < Ti < T c , the burning ceases and the site becomes passive again. Part of the heat release is 
lost (radiated away), while the other part, aJ, is transferred (isotropically) to the closest neighbors. 
Summarizing the above, the heat release can be written approximately as 



where 8{x) is the step function. The term 8(—dT/dt) is an approximate manifestation of the his- 
tory dependent (hysteresis) burning for T c > T > T\ (burning now if it was burning at the previous 
moment/step and not burning otherwise). This expression is not quite correct for the temperature 
of a site does not have to change monotonically when an avalanche develops. We leave the more 
detailed discussion of this for another paper especially devoted to this model. For the purposes 
of the present discussion such details are irrelevant, and we consider (ITOb as a sufficiently precise 
description of the burning process. If an active site would be left alone, its temperature would de- 
crease as T = T(0) exp(— fit). Here the quantity r k (1//j) m(T c /TJ) has the meaning of the life 
time of an active site if it were not affected neither by other sites nor external driving. Let A t be 
the time step and A/ the site size. The amount of heat a site x is receives is given by 



dT 
~dl 



qp(l-r)(t)), 



(9) 



J = fiT9(T c - T)6(-dT/dt)8(T - TJ) + fiT8(T - T c ) 



(10) 



dT(x) 
dt 



qp[l - v(t)} + a[J(x + Aj) + J(x - Aj)] - J(x), 



(11) 



which we write in the following form 



dT(x) 
dt 



qp[l-r]{t)} + (2o- 1)J + 



oAj d 2 J 
2 dx 2 ' 



(12) 
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Integrating (fTTT) over a cluster of the size w, one gets 

— J Tdx = qpw + (2a — 1) J Jdx — J&, (13) 

where we averaged over time the random fluctuations of the input r\. The last term is the heat flux 
at the boundaries. 

The probability of growth should be proportional to the heat flux from the active site at the clus- 
ter boundary to the neighboring passive sites. This probability should depend on the temperature 
of the passive sites. In the stationary regime the time-average growth probability would be de- 
termined by the average temperature T p of passive sites. Thus, growth is essentially independent 
of the cluster size. Respectively, the shrinking probability depends on the state of the boundary 
site and is not particularly sensitive to the cluster size either. In this case the parameters a and (3 
are constant, and one expects that the cluster distribution is an exponential, / oc exp(— w/wq). 
However, if the heat transfer in the active area is suppressed (active sites do not easily accept heat 
from active neighbors) spreading from the central regions with the constant speed up to the cluster 
boundaries, one estimates that P + oc l/w, while P_ w const. In this case a = A(ai — a 2 /w), 
and (3 = (ai + a 2 /w)A 2 /2, and / oc (w + w ) x exp(—w/w c ), where w Q , w c and A are constants. 
In the range w <^ w <^ w c (if such this range does exists at all) a power-law distribution should 
be observed. In the opposite case, when the heat is transferred immediately from the inside to the 
cluster boundaries, P + oc w and P_ rs const, no stationary state can exist, since a < for suffi- 
ciently large w. Such systems are unstable and are disrupted into avalanches which will cover the 
entire system. 



IV. CONCLUSIONS 



We proposed a kinetic approach to the description of avalanching systems, defining a distri- 
bution function f(w,t) for the clusters of active sites. In this way we derived a kinetic equation 
for the temporal evolution of f(w, t) and analyzed its steady state limit in the inertial range, suffi- 
ciently far from the smallest scales where driving explicitly shows up, and sufficiently far from the 
larges scales where fractality and merging become progressively more important. The stationary 
distribution function f(w) depends, in general, on the probability of the micro-processes resulting 
in cluster growth and shrinking, that is, the processes governing the evolution of avalanches. In the 
case of locally induced growth at the boundaries the shape of the distribution is determined by the 
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dimension of the system (or fractal dimensions of clusters if they are not compact). There is no 
sensitivity to the input details. The obtained universal shape of the distributions is not limited to 
the weak driving regime or to the system dimension above some critical value, and can be used for 
direct and easy comparison with experiments and numerical modelling. The total average driving 
should affect the state of the system, as we have shown in a particular model. The estimates given 
in the present model represent just the first step toward a more elaborated kinetic model of the 
dynamics of avalanches. We remark that our analytical predictions have been checked by ID and 
2D burning model simulations to be reported elsewhere. 
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